New approximation for nonlinear evolution in periodic potentials 
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A new approximation for evolution described by Nonlinear Schrodinger Equation (NLS) with pe- 
riodic potential is presented. It relies on restricting dynamics to one band of the bandgap spectrum, 
and taking into account only one, dominating Fourier component in the nonlinear Bloch-wave mix- 
ing. The resulting equation has a simple, discrete form in the basis of linear Wannier functions, and 
turns out to be very accurate as long as the modes in other bands are not excited and the potential 
is not very deep. Widely used approximations, the tight-binding approximation and the effective 
mass approximation, are derived from the new equation as the limiting cases. 
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The Nonlinear Schrodinger Equation (NLS) is one of 
the most common and basic nonlinear wave equations. 
Recent experimental progress in the fields of nonlinear 
optics and Bosc-Einstein condensation, where it is known 
as the Gross-Pitaevskii equation, has strongly stimu- 
lated theoretical studies on its properties and solutions. 
The ability to describe and understand nonlinear evolu- 
tion in periodic potentials is of a great practical impor- 
tance, since this kind of potential can be easily created 
by a periodically varying refractive index of an optical 
medium [H, H, d, 0|i or by a standing laser wave im- 
posed on a Bose-Einstein condensate @, Q ■ On the other 
hand, the presence of a periodic potential can dramati- 
cally modify the wave diffraction, and creates a way to 
control the nonlinear evolution of physical systems. 

To date, there are a few models commonly used to 
describe these systems in a simplified way. In the case 
of a deep potential, the essential properties of system 
dynamics are included in the tight-binding (or coupled- 
mode) approximation, which assumes that the interac- 
tion of neighbouring cells can be described by tunneling 
between their localized modes. As a result of this as- 
sumption, one obtains a Discrete Nonlinear Schrodinger 
Equation (DNLS), a very convenient tool for studying 
and modelling e.g. light propagation in weakly coupled 
waveguide arrays [H, 0, 0] • If the Bloch-wave spectrum 
of the system contains only a small part of the Brillouin 
zone, one can use the effective mass approximation to de- 
scribe motion of a wavepacket spanning over many sites 
of the potential It was successfuly used to model 
the creation of gap solitons - nonlinear stationary states 
exisiting even in the case of a negative (self-defocusing, 
repulsive) nonlinearity 0, @] • 

In this letter, a new approximation is presented, which 
brings together positive aspects of these two models, 
while preserving most of their simplicity. The derivation 
relies on assuming that the dynamics restrict to one band, 
and taking into account only one, dominating Fourier 
component in the nonlinear Bloch-wave mixing. It leads 
to a simple discrete equation in the basis of linear Wan- 
nier functions 8, 9]. The new equation can be used to 
describe evolution in weak or deep potentials (as long as 



they are not very deep), and includes dynamics in the 
whole Brillouin zone. Its accuracy is better than could 
be expected, as will follow from direct numerical compar- 
ison of its results to solutions of the NLS equation. 

In reduced units, the Nonlinear Schrodinger Equa- 
tion with periodic potential for time-dependent function 
u(x, t) reads 



,du d 2 u , . l2 



(i) 



where U(x + 1) = U(x) E M and a = ±1 is the 
sign of the nonlinearity. The norm of the solution 
A" = \u\ 2 dx is a constant of motion. After per- 

forming Fourier transform one obtains the equation for 

u(k, t) = (2tt)- 1 /2 u ( a . j t ) e - ikx dx 



.du{k) 



= -k 2 u{k) - U n u(k - 2im) 



(2) 
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u(k\)u* {k2)u{k — k\ + fc2)dfcidfc2 , 



where U(x) = J2 n U n e 2l ^ nx , and the integer index n £ 
(— oo, +oo). 

From now on the dynamics will be restricted to one 
certain band of the bandgap spectrum. The pseudomo- 
mentum will be denoted with k' . For each k' € (— 7r, tt), 
there exists a Bloch function v(x, k') of energy E(k'), and 
its Fourier transform v(k, k') is the solution of the linear 
version of Eq. (f2]) 

E{k')v{k, k') = k 2 v(k, k') + Y u nv{k - 27m, k'), (3) 

n 

with orthonormality condition J v(k,k' 1 )v*(k,k' 2 )dk = 
S(k[ — k' 2 ). It is assumed that the solution is the su- 
perposition of Bloch functions of the chosen band 



u(k,t) = / (Hfc'^MMOdfc' 



(4) 



This assumption implies that the nonlinearity is rather 
weak, since strong nonlinearity leads to Bloch wave mix- 
ing of waves in different bands, or even destruction of 
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FIG. 1: (color online) a) Linear Bloch wave spectrum in 
cosine potential of depth e = 4. Small dots correspond to 
frames in Fig. [2] b) Contribution of the main component 
in decomposition of a Bloch wave in Fourier space. Dashed 
and solid lines correspond to the first and the second band, 
respectively, c), d) Same as a), b), but for potential depth 

£ = 10. 



the band structure [5j. After substituting this equation 
into ([2]) and using ([3]), one can multiply both sides by 
u*(k, k'o) and integrate over k to obtain 

-f-oc 

v*(k^ 1 k' )v(k ll k' l )v*(k 2: k 2 ) x 

■oo 

x i/(kQ — ki + k 2l k' 3 )dkodkidk2- 

For this equation to be useful, the nonlinear term has 
to be simplified. According to Bloch theorem, Fourier 
spectrum of a Bloch function is discrete, u(k, k') — 
J2n l/ n(k')5(k — k' — 27m). It turns out that in most 
cases one of the Fourier components dominates, as can 
be seen in Fig. [TJ where the results for cosine poten- 
tial U (x) = — e cos(27ra;) are shown. Let riband be the 
band number, starting from zero for the lowest one; then 
one can assume that v{k,k') w 5{k — fcmain)) where 
&main = k' + signum(fc')7m ba nd for even n ban d, and 
femain = k' - signum(/c')7r(n ba nd + 1) for odd n ban d- This 
property is related to the fact that in the limit of van- 
ishing potential Bloch decomposition is equivalent to the 
Fourier transform. From Fig. Q] one can see that this as- 
sumption is fulfilled better for a weaker potential, and 
that the agreement is worst at the edges of the Brillouin 
zone and in its center for higher bands (these are the 
points of jump in fc ma i n )- However, as it will be shown 
below, this approximation gives good results even in the 



case of gap solitons, nonlinear objects which reside in the 
vicinity of these points. 

As a result, the above equation takes the form 



1 ^ = -^')#' ) i) + 



dt 



2n 



<f>(k[,t)<f>*{k 2 ,t)cf>[k' -k[ + k' 2l t)dk[dk 2 . 



Inverse Fourier transform gives a simple discrete equa- 
tion, which is the main result of this letter 



dt ^ 



(5) 



Here cj>(k',t) = ^tt)" 1 / 2 £ (j)„(t)e- ik ' n and E(k' 



Ik' I 



, E- n = E*. The norm of the discrete func- 



tion is "Y^ n |</>n| 2 = N. Most of the terms in the infinite 
sum can be neglected, since absolute values of E m are 
significant only for m close to zero. In numerical sim- 
ulations presented below, only terms with \m\ < 5 were 
taken into account. Values of E m can be found quite eas- 
ily, by solving a linear eigenvalue problem (|3|) (see e.g. [5() 
and performing an inverse Fourier transform. 

The scheme of the above derivation is similar to the 
one presented in [8j , where a vector discrete equation was 
obtained for evolution in the basis of Wannier functions. 
In fact, the one-band assumption imply that the function 
u(x,t) can be approximated by 



u(x 7 1) w u$(x, t) = 2J <t>n{t) w ( x - n ) 



(6) 



where w(x) is the linear Wannier function w(x) — 
f* n i/(x,k')dk', see Fig. |2^l). Therefore, <p n can be in- 
terpreted as the amplitude of the wavefunction in the 
n-th site of the periodic potential. 

It is relevant to point out possible generalizations of 
the presented model. The one-band approximation can 
be extended to include more than one band. In this case, 
one would obtain a system of discrete equations , each 
of them coupled with others by cross-phase modulation 
terms. The derivation can also be easily generalized to 
the case of multi-dimensional NLS equation 0, [HI • How- 
ever, one has to be aware that in this case the nonlinear 
mixing between bands is more likely to occur, since the 



band gaps are not always closed Finally, the pre- 

sented method can be applied to other types of nonlincar- 
ity, e. g. quadratic, cubic-quintic, or other nonlinearitics 
after expanding them in Taylor series. A detailed study 
on these generalizations will be presented elsewhere. 

Equation ([5]) has two interesting limiting cases. For 
a deep potential, the energy dependence E(k') becomes 
close to the cosinus function, and the infinite sum can 
be approximated by three terms with m = — 1, 0, 1. This 
form of the equation is equivalent to the one obtained 
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FIG. 2: (color online) a) Localized Wannier function of the 
first band (thick line) and the second band (thin line) for 
cosine potential depth e = 4. Dashed line shows the shape 
of the potential, b) Stationary soliton state at the bottom of 
the first band (see Fig. []J,) , according to the full NLS model 
(Eq. [T] solid line) , and the approximate model (Eq. [5] dashed 
line, the two lines overlap). Dots show the amplitude of the 
discrete function <f> n . Parameters are: e = 4, a = 1, N = 3.4 
and j3 = — 1. c) Soliton at the top of the first band for a = — 1, 
N — 7 and (5 = 8.5. d) Soliton at the bottom of the second 
band for a = 1, TV = 6.9 and (3 = 11. 



within the tight-binding approximation [J, 2, HI . How- 
ever, the new equation will not be accurate, since for deep 
potentials the assumption made in the derivation is not 
well satisfied, see Fig.[T]i). Still, it will give qualitatively 
correct results. 

Another limit corresponds to the case when the 
function (f> can be written in the form <j) n {t) — 
<&(n,t) exp(ik' Q n), where $(x,t) is an envelope slowly 
varying in the spatial coordinate x. Using this relation 
and expanding <f>(x,t) in Taylor series up to the second 
order leads to the effective mass equation 

i _ = _ EW)4 _ j „ ! _ + __ +H ^ 4 , (7) 



where v g = (dE/dk')\k' 
m cS = (d 2 E/dk' 2 )- l \ k , =l 



■k' is the group velocity, and 
' is the effective mass. In this 



case, the new equation (|5]) agrees quantitatively with the 
NLS equation |T]), This is confirmed by numercial simu- 
lations, see Fig. O 

To test the new approximation, the equation J5|) 
was applied to description of band-gap solitons [3, Q. 
These states are the inherently nonlinear solutions of the 
NLS equation fl} in the form u(x,t) = u(x) exp(— ifit), 
where (3 is the eigenvalue of the soliton, lying in 
the gap of the band-gap spectrum. In Fig. [2] these 
states are compared with analogous solutions of Eq. (J5]), 
4> n {t) = 0„exp(— if3t) and the corresponding functions 
it0, cf. Eq. In particular, the two figures a), b) 




FIG. 3: (color online) Comparison of the calculated soli- 
ton width (dark lines) and maximum amplitude (light lines) 
according to the NLS equation ([T]) (solid) and the approxima- 
tion ([6]) (dashed), versus the potential depth e. Structure of 
solitons is the same as in in Fig. (2b). For each e, the norm of 
the NLS solution is TV = 3, and the eigenvalues j3 of the two 
compared solutions are equal. The vertical line indicates the 
value of e when the widths of the first band and the first gap 
are the same. This is the point of transition from the weak 
potential regime to the deep potential regime. 



present solutions of the same equation, describing evolu- 
tion in the lowest band, with the only change in the sign 
of the nonlinearity. The agreement with the full NLS so- 
lution is perfect. Interestingly, the solution in Fig. ^p) 
has larger width and norm than the solution in Fig. [^)- 
In the tight-binding model, solitons with larger norm al- 
ways has a smaller width, independently on the sign of 
nonlinearity [f3 |. Here, the new equation takes into ac- 
count the difference in diffraction strength at the top and 
the bottom of the first band. In Fig. a soliton com- 
posed of Bloch waves from the second band is shown. In 
this case, maxima of amplitude lie on maxima of the pe- 
riodic potential. The agreement between the full model 
and the approximate model is in this case somewhat infe- 
rior. Additional simulations have shown that the reason 
for this is the strength of the nonlinearity; the eigenvalue 
of the solution is shifted deeper into the band gap than 
in the previous cases. In general, the equation (0 works 
best if the nonlinear energy is much smaller than both 
the gap width and the band width. 

In Fig. [3] a systematic comparison of the full and ap- 
proximate models is presented. Here the width and max- 
imum amplitude of the lowest-energy soliton is depicted 
versus the potential depth e. Good agreement of cal- 
culated width W(u) = 3/ |a;||M(x)| 2 dx / J \u(x)\ 2 dx sug- 
gests that the shape of the approximate solution is correct 
even for very strong potentials, with the only difference 
in the norm N. 

The new approximation can also desribe dynamical 
evolution problems. A very simplified version of it was 
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FIG. 4: (color online) Motion of a soliton wavepacket with 
imprinted linear phase for N = 3 and e = 4, according to 
a) the NLS equation and b) the discrete equation c), d) 
Same as in a), b), but for e — 10. In this case, the soliton 
is trapped by the Peierls-Nabarro potential. The shape of 
solitons is the same as in Fig. [2p). 

already used for justification of the phenomenon of spon- 
taneous migration of Bloch waves to the regions of normal 
diffraction for positive nonlinearity, and to the regions of 
anomalous diffraction for negative one |7] . This nolinear 
phenomenon was observed experimentally in two differ- 
ent systems Il3| . and can be utilized for efficient soliton 
generation [7j. Here the new model is used to describe 
soliton mobility, see Fig. [4] The lowest-energy soliton 
wavepacket was "boosted" by imprinting a linear phase 
u(x,0) = u{x) exp(ikx), with k = 0.1. fn the case of 
a weak potential, the wavepacket started to move with 
a constant velocity, fn the case of a deep potential, it 
has been trapped in central sites by the Peierls-Nabarro 
potential [l4j. This effect would not be seen within the 
usual effective mass approximation. Here, the agreement 
between the NLS equation and Eq. ([5]) is very good for 
the weak potential, worse for the deep potential, but the 
main effect is still apparent. 

In conclusion, a new approximation for evolution de- 
scribed by Nonlinear Schrodinger Equation with periodic 
potential was presented. The derivation is based on the 



one-band approximation and simplification of nonlinear 
Bloch-wave mixing, and leads to a simple discrete equa- 
tion in the basis of linear Wannier functions. The equa- 
tion works very good as long as the nonlinearity does 
not cause excitation of modes in other bands and the po- 
tential is not very deep. The new model was used for 
description of gap solitons and a dynamical evolution. 
It was shown that the tight-binding approximation and 
the effective mass approximation can be derived from the 
new equation as the limiting cases. Possible generaliza- 
tions were pointed out. 
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